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Abstract 

The problem of construction of the surfaces with given sets of the 
points with horizontal tangential planes is considered. Such consider- 
ations are of interest in the problem of computer simulations of the 
waved ocean surfaces. 



1 Statement of the problem 

1.1 Physical statement 

Various problems of optical/acoustical monitoring of wavy water surface are 
related with simulations of the surfaces with given properties [TH5]. The 
example of it is the problem of construction of two-dimension surface having 
horizontal target planes at the given points (Luchinin A.G., 2003). Variants 
of solutions of this problem are presented in this article. 

1.2 Mathematical formulation 

The main problem A. Let DJt = |xj, i = 0, V} be the finite subset of the 
given open set Q C -R",n > 1. Construct a function F (x) of the given class 
^ & (Q) possessing zero gradient VF (x) = at all points of VJt . 
We study also the next specifications of the main problem: 

B. the values of the function at the points DJt are given (F {'Xi) = Zi, i = 
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C. gradient of the function takes zero values only in DJl. 

It is clear that sets of the solutions of A-C are infinite. In general, 
solutions of the main problem A do not solve problems B,C. The main goal 
of the next research is to determine the existence of the solution of a problem 
in a given class of functions ^. 

2 Analytical approach 

Let ^ = C^{Q), 1 < k < oo. For this class of functions we can easily 
construct an analytic solution of the main problem. One of possible con- 
structions is described below. Let |-| be Euclidean norm in R^. Consider for 
every point Xj G OJl the open sphere 

S'j = {x : |x — Xj| < Vi} C Q 

with the radius rj > and center point Xj. Let > 1, q be arbitrary values, 
(fi (■) G C'^ (-R") be arbitrary functions such that cpi (0) = 1, \/(pi (0) = 0. 
Let 



Xi{^)= < fli ~ \ |x| < Ti, 
[ 0, |x| > Ti, 

N 

^ (^) = (X - X,) (pi (X - Xi) . (1) 

i=0 

It easy to prove the following statements. 

Proposition 1. For every set of sufficiently small values ri (such that 
Xj ^ Sj, i ^ i) and for any values Ci the function (QP zs a solution of the 
problem A. In addition, if 

1. Ci = Zi, i = 0, A^, then (CP gives a solution of the problem B; 

2. all spheres Sj do not intersect, then ^\) gives a solution of the problem 
C. 

Proposition 2. Let n = 1, xq < Xi < . . . < x^ . Then under the 
conditions 

Ti < \xi - , Ci^iCi < 0, i = 1,N, 



2 




Figure 1: The surface ([T]) for the set of 10 random points marked on the 
graph by stars. 

the function (QP gives a solution of the problem C. 

The pecuharity of the described construction is the provided by multi- 
phers Xi (x) finiteness of the terms of the sum ([T]) ; this provides zero values 
of gradient at the points of 971. Note that using the corresponding values of 
gradients V^pi (0) we can assume these gradients not only zero but also any 
given, so we can use the construction ([1]) for solving of more complicated 
problems to construct of surfaces with given slopes at the given points. 

From physical point of view the solutions ([T]) can simulate only oscillation 
of too viscous liquid. An example of the numerical simulation by means 
of formula ([1]) for = e, (pi (x) = 1, 10 random points 971 and random 
coefficients Cj, is presented on Fig. [H 

3 Spline approach 

In this section we use splines of various types to solve problems A-C for 
n = 1, 2. Polynomial splines of the 4-th order (see sec. I3.1.H I3.1.2p make it 
possible to solve problems A and B. As for the problem C then in the case 
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n = 2 spline approach is useless if only the set DJl does not have the form of 
the product 

But in the case n = 1 the solution of the problem C can be found in the 
class of trigonometric sphnes (sec. I3.2.ip . 

3.1 Polynomial splines 

3.1.1 One-dimensional problem (n = 1) 

Consider the following problem B: for given sets a = xq < Xi < . . . < Xn = b 
and {zi,i = 0, A^, } build a spline of the 4-th order F (•) G [a, b] such that 

F (x) = Si{x) , X e [xi-i,Xi] , i = 1,N, 

Si (x) = ai + bi {x — Xi) + Ci{x — Xj)^ /2\ + di{x — Xi)^ /3! + Cj (a; — Xi)'^ /4!, 
F (xi-i) = Zi^i, F {xi) = Zi, 
F'(xi_i) = 0, F'{xi) = 0. 

(3) 

Now we prove that such splines exist. We have for the coefficients of the 
spline 

(a^i) Q-j Zi, 
F' (xi) = bi = 0, 

< F = ai + biAxi + CiAx]/2 + + eiAxf/24 = Zi_i, 

F' = bi + CiAxi + diAx'i/2 + aAx^Q = 0, 

^ Axi = Xi^i - Xi, i = 1, N. 

So we obtain the equations 

Ci = —diAxi/2 — CiAxf/G, 
-diAxf/l2 - eiAxf/24: = Azi = Zi^i - Zi, 
di + CiAxi/S = —2ci/Axi, 
di + eiAxi/2 = -12Azi/Axf, 

and we can evaluate all coefficient by q, 

r a, = z„ 

bi = 0, 

< Ci = Ci, (4) 

di = -6 (ci - 4:AziAx~'^) /Axi, 
^ ei = l2 {a - GAziAx;'') /Axi 
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Spline and its derivative are continuous both. In general, for arbitrary co- 
efficients Cj we obtain smooth splines with discontinuous second derivative . 
The examples of such splines are presented on Fig. |2] 
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Figure 2: Examples of smooth polynomial splines of the 4-th order solving 
the problem B. The choice rule for the coefficients : Cj = ^, where ^ is 
a random value uniformly distributed in the interval [—0.05,0.05]. Points 
i) are marked on the graphs by circles. 

If we do not fix parameters Zi in (jl]) then splines solve only the problem 
A (see Fig . [3D 

C^-spline Let us now construct a spline with the continuous second 
derivative. For realization of it we must require 

<_i(a;,_i) = <(x,_i), i=TJI, 
that tends to the equations 

Cj_i = Ci + diAxi + -^^x'^i = Ci- 12AziAx~'^. 
This provides the recurrence 

a = Q_i + UAziAx-^, (5) 
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so 



Q = Co + U^AzjAxJ^, i = l,N. 

Hence we obtain a one-parameter aggregate of splines of the class [a, b] . 
The parameter of it is cq, the second derivative of the spline at the initial 
point a. By choosing of cq, we can construct splines with different properties. 
For example, when 



-2 
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Co = —12 min "^^^AzjAx^ 
all points of 9Jt are points of local minima (see Fig. HI a), for 



Co = — 12 max ''^^AzjAxj ^, 



all points 9Jt are points of local maxima (see Fig. HJ b), and for 

12 . ^ ^ 
Co = y'AzjAx'^. 

3 = 1 
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Figure 4: Polynomial splines of the forth order with additional properties: 
a) all points of DJl are points of local minima, b) all points of DJi are points of 
local maxima, c) the mean value of the second derivative for over all points 
of dJl takes its minimum. 

the mean value of the second derivative for over all points of DJt takes its 
minimum (see Fig. HJ c) 



3.1.2 Two-dimensional problem (n = 2) 

Consider the following two-dimensional problem B: for given set of points 
{xk, yi, Zki), a = xo < xi < . . . < xk = b, c = yo < yi < . . . < yi = d build a 
spline of the 4-th order F (■) G ([a, b] x [c, d]) such that 



Ski [x, y) = Y1 ^^iT' - {y - yiY 



F (x, y) = Ski {x, y) , (x, y) G Pm = [xk-i, Xk] x [yi-i,yi] k = 1,K, / = 1, L, 

4 

{x - XkY {y - yiY 

i,j=Q 

F {xk,yi) = Zki, 
[ VF{xk,yi) = 0. 

(6) 
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By decomposition of the function into one-dimensional splines onto 
grid lines {xk, yi) ,k = 1, K, I = 1, L, we receive the following equations for 
the coefficients of the two-dimensional spline, 



i=0 \j=0 
4 



i \{x- Xkf = 

i=0 



(k l-l) 



x-Xk) , y = yi-i, 



y = Vh 



i=0 



i=0 



X^/3f ^^Hy- yiY , X = Xfc„i, 



(7) 



j=0 \i=0 
4 

Xl^of (y - yi 

3=0 



j=0 



yiY, 



X Xk^ 



where we know values 7^ , /3 j from formulae (111); so we know the first col- 



umn and the first row of the matrix A^^^^ = I a^f^ \ . Last 16 coefficients 

I ^ ) i,i=0 



zu /3f') 

{kl) (kl) (kl) (kl) 

'^ll ^12 ^^13 ^^14 

{kl) (kl) (kl) (kl) (kl) 

72 '^21 '^22 '^23 '^24 

(kl) (kl) (kl) (kl) (kl) 

^3 ^^31 ^32 ^^33 ^^34 

(kl) (kl) (kl) (kl) (kl) 

74 ^il '^42 '^43 '^44 



are unknown. For them we have from ([7]) the equations 



(fc-1 i) n{kl) 



i 5 ^ -'-5^) 



4 



(8) 



A/3f^ J = 1,4. 



v> i=i 
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The matrix of this system is 
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and its rank is equal 7. If we choose 

{kl) {kl) {kl) {kl) {kl) {kl) {kl) 
^11 1 ^12 1 ^13 1 ^14 1 ^21 1 ^31 1 ^41 1 

as known values (they correspond to nonzero minor of the order 7) and choose 
other 9 elements as free, then we obtain from ([8]) the system 

' a!f Ay, + a[fAyf + a[fAyf + a[fAyf = A^['^^ = 0, 
a^f Ay, = A7f - Ayf + a(f Ayf + a^f Ayf , 
4f Ay, = A7f ) - 4? Ayf + 4f Ayf + 4f Ayf , 
< 4f Ay, = A7f ) - 4? Ayf + 4f Ayf + 4f Ayf , (9) 
aifAxk + aifAxl + af^xl + 4f Aa;^ = A/3f ^ = 0, 
a[fAx, = A/3f - a^fAx^ - 4f Ax^ - 4f Ax^ 

^ a[fAx, = A/3f - ag'^Axl - ag^^Axl - a^^ Axl 
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with the solution 



a 



(kl) 

13 
(kl) 



.(kl) ^ 



2 ^^k' 



' = Ap. 

jikl) _ ^Jkl) 

4f = A.r^Ay-^ 



(fcO A 2 _ ^C^O A 3 
"33 ^"''fc "43 ^-^fc' 



„v ; A 2 _ ^C^O A 3 
"32 '-^•^k "42 ^"''fe' 

"i?^2/f - a[fAyf, 



— a 



"21 ' = ^7?' ^ At/; 

(fcZ) (fc/) 

"ll = ~" 



(kl) 

42 
(fcO 
^32 



a 



(fcZ) 
14 



21 

(kl) 



Axk - 



a\r>Ay 



- "22 - 

"if Ax^ - 

- "SfAy,- 



(H) 
^33 

(kl) 
^23 



Ayf 
Ayf 



(H) 
^34 

^24 



Ay\ 
Ay\ 



(kl) 

'■13 



Ayr'- 



(10) 



C'^-spline Let us set, for example, in ([101 



2,4,j = 2,4. 



Then we obtain from flTU]) 



^13 

(W) 
''12 

(W) 
Ml 

^31 
(kl) 

H\ 

(kl) 
n\ 
(kl) 

*14 



A/3; 



A/3f')Ax,, 

A|/z, 

Air^Vu 
Ai?^Am, 



'2 

A^r 

(kl) , 



{kl) 



A7. 
"If Ay, 



+ A7f)Axfc + A74 



A/3f)Ax,^Ayr^ 



('^'^Ax? ] AxfcA?/,, 
A/3j 



^(ki) 



(11) 



The spline with such coefficients is continuous but, in general, its first deriva- 
tives have jumps on the boundaries of the cells Pfc^. The example of such 
spline see on Fig. |5] 

C^-spline Let us try to choose such free parameters that spline ([6]) 
becomes C^-function. It is sufficiently to require that the first derivative 
dF/dx is continuous at all vertical lines x = x^ and the dF/dy is continuous 
at horizontal lines y = yi of the grid. This tends to the equations 



(kl) 



E 



ta. 



1=1 



(fc+i 



As^fc+i? "ii 



{kl) 



2, 4, A; 



l,is:- 1, / = 1,L- 1. 



(12) 
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Figure 5: The example of polynomial C°-spline with the coefficients ( ITTl) . 
The points of Wl are marked by stars, K = L = 10. 



We see that the second row and the second column of the matrix are com- 
pletely determined by the coefficients of the matrix from contiguous cells 
(that lie to right and above from a given cell). On the other hand due to the 
continuity of the spline and the equations ffTOl) these elements can be found 
by nine coefficients of the matrix 



"22 
"32 
"42 



a 



(kl) (kl) 



23 



24 



a 



(kl) (kl) 



33 



34 



a 



(kl) ^{kl) 
43 



'44 



(13) 



From the first and the fifth equations of the ([9]) and from the (fT2|) we receive 



4 4 



0, 



4 



j=l i=l 
4 4 



J a. 



i=l 



i=l j=l 
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So we obtain for a- 



(ki) 



/ 4 4 



1=1 j=l 

4 4 



0. 



j=l i=l 

Five additional equations we obtain from the system 



4f + aifAx, + alf Ax^ = A^t^Axf - a^f Ax,-\ 
ag) + a^f Ay, + a^f Ayf = A^^^Ayf' - a^f Ayr\ 
+ aifAy, + aifAyf = A^f^Ayf - Ay-\ 
ai? + aif A^, + alf Ayf = A^f ^ A|/r^ - a'^ Ay-\ 



A^r^Ay-^ - 



(14) 



(15) 



From systems f|T5|) and f|T4|) we receive the system of seven equations for nine 
coefficients (fT3l) 



a 



(feO 

22 



^ "32 



Axfc + a^io'^Ax^ 



A42 



^23^'* + ai^s^Axk + a)^^'Ax 
"22 



a 



(ki) 

32 

(feO 

^42 



^ "23 



"33 

-t- 



Ay, + a^f Ayf 
Ay, + aif Ay2 

44 '-^iJi 



EE 

i=2 j=2 
4 4 

EE 

I, i=2i=2 



j4" 



^yi + a 
Ayi-'Axt' 



talf^Ayj '^Ax 



k 



- A(3i''^Axf 

- A^f^Axf 
A^'^Ayl- - 

M^^vi' - 

A7f)A.,-2 

4 



(ki) 

''12 



Ax 



k ' 

"13 ^^fc ' 



(M) 
*21 

(H) 
'-31 

(fc/) 
Ml 



1 



Ay, 



-Y^ja^^Ayr - J:-^i'Ax, 



i-l 



3 = 1 
4 



i=2 

^^ai^^Axt\ 



i=2 



(16) 



12 



However the matrix of this system 



1 





1 








2AyiAxk 


2AxkAyi 


Axk 








1 





2AyiAxl 


SAxlAyi 
4 Axl Ay I 


Axl 











1 


2AyiAxl 





1 


^yi 








3AyfAxk 


2AxkAyf 
SAxlAyf 
4Ax|Ay? 
2Aa;fcAyf 
3Aa;^Ayf 





Axk 










3AyfAxl 





Axl 








^yi 


SAyfAxl 








Ay? 








AAyfAxk 











Ay? 





AAyfAxl 














Ay? 


AAyfAxl 


4Ax|Ayf 



has the rank 5. This means that system f|T6l) has no solution. Hence two- 
dimensional problem B has no solution in the class of C^-splines. 

Besides the problem of smoothness the next problem is the appearance of 
false stationary points that do not belong to The above considerations 
show that class of polynomial splines is not convenient for solutions of the 
problem B,C. 

3.2 Trigonometric splines 
3.2.1 One-dimensional problem 

Consider the following problem B: for given sets {xi,Zi), a = xq < Xi < 
. . . < Xn = b build a trigonometric spline F {■) E [a,b] such that 



Conditions (fT7|) do not identify the spline. For example, we can take the 
following coefficients 



tti = {zi^i + Zi) /2, 
h = {zi_i - Zi) /2, 

< Q = Tt/ A Xi, 

di = -iTXi-i/ A Xi, 




F (x) =ti{x) , X e [xi-i, Xi], i = 1, N, 
ti {x) = ai + bi cos (cjX + di) , 

F {Xi) Zi, 



(17) 
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The corresponding spline has the form 

. / \ ^i— 1 , ■^j— 1 
ti[X) = COS TT . 

^ ^ 2 2 Ax, 

For such coefficients false stationary points are absent; this means that really 
the spline solves the problem C. 

In general, for arbitrary values the spline has jumps at points Xj. If to 
forbid the jumps, 

t'U i^i-i + 0) = -vr^ {z,_2 - Zi-i) A x-y2 = 
= TT^ - Zi) A xry2 = t'l {xi-i - 0) , 

then we receive the relations 



{zi - Zi^ / A = ;U (-1)* , 2 = 1, n, 

where /i is an arbitrary value. Such kind ^-splines form two-parametrical set 
as for fixed value zq other values z, can be calculated by the formula 

i 

z, = z, + ^lY,^-l)' /\xf (18) 
1=1 

which contains two free parameters zq and ^. For fixed set {xj} the parameter 
/i evaluates mean deviation of the spline from the initial value zq: 

n n 

J2 - ^o) /n = fiJ2 (-1)' (1 - ^ 

i=l 1=1 

Examples of trigonometric splines (fT7|) are presented on Fig. [6] 
3.2.2 Two-dimensional problem 

Consider the following problem B: for given sets {xk,yi, Zki) , a = Xq < 
xi < . . . < xk = b, c = Ho < Hi < . . . < i/L = d build a spline F {■) e 
{[a,b] X [c,d]) such that 



F {x, y) = tki {x, y) , (x, y) G P^, k = l,K, I = 1,L, 

hi {x, y) = [ttki + bki cos (cfe/x + 4/)] [Aki + Bki cos {CkW + Du)] , 

F {xk,yi) = Zki, 

yF{xk,yi)=0. 
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Figure 6: A beam of smooth trigonometric splines. The values {zi} were 
taken randomly from a given interval. 

In analogues to one- dimensional case the conditions (IT^ do not identify 
the spline. For example, we can take 

Cki = vr/Axfe, 

dki = -TTXk-i/Axk, ^20) 
Cki = n/Ayi, 
Dki = -rcyi_i/Ayi, 

where Axk = Xk - Xk-i, Ayi = yi - yi-i. In this case S/F{xk,yi) = at 
the vertex of Pm (and only at these points; so, actually the spline solves the 
problem C). For other coefficients of the spline we have the equations 

tki {xk-i,yi-i) = [aki + bki] [Aki + Bki] = Zk-u^i, 
tki {xk,yi-i) = [ciki — bki] l-Aki + Bki] = Zki-i, 
tki {xk-i, yi) = [ciki + bki\ [Aki — Bki] = Zk-iu 
tki {xk, yi) = [aki - bki] [Aki - Bki] = Zki 
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which give 



dkl^kl 
bkl^kl 
0-klBkl 

bkiBki 



{zk-ii-i + Zk-ii + Zki-1 + Zki) /4, 

{Zk-ll-l + Zk-ll — Zkl-i — Zkl) /4, 

{zk-ii~i + Zki-i — Zk-ii — Zkl) /4, 
{zk^ii^i + Zkl - Zk-ii - Zki-i) /4, 



(21) 



that is sufficient for the sphne determination. As coefficients of tki {x, y) 
depend only on the sphne values at vertexes of P^;, then we must restrict 
these values in order to provide the smoothness of the spline. 

It is easy to prove that we need to provide the smoothness of the spline 
on the boundaries of Pm only. Take for example the common boundary of 
Pki and Pfcz-i, the segment [xk-i,Xk] x {yi_i}. The spline is continuous on 
this segment if and only if for all x G [xk~i,Xk] 



Substitute formulae (I2T1) to these equations we receive the identities for 
Zkl- Also we get the identities when check the continuity of the spline on the 
vertical boundary {xk-i} x [yi^i,yi]. Hence the spline is continuous for all 
values {zki}. Actually the spline belongs even to the class {[a,b] x [c, o?]). 
To prove it we must prove the continuity of the gradient on the boundaries 
of Pki- For example, using fl22]) . we have on the segment [xk-i,Xk] x {yi-i} 

F^(x,i/;_i±0) =0, 

{x, yi-i ± 0) = nbki sin (tt (x - Xk-i) / Axk) [Am + Bm] / Axk = 
= Tibki-i sin (tt (x - Xk-i) / Ax^) [Am-i - Bki-i] / Ax^. 

If we want to construct the more smooth spline then we must to impose its 
values {zki}- For example the condition of continuity of the second derivative 
F'a on the boundary [xk-i,Xk] x {yi-i} tends to the equations 



[au + bki cosTT (x - Xfe_i) /Axfc] [Au + Bki] = 
[ttki-i + bki^i cosTT (x - Xfc-i) /Axfc] [Aki-i - Bki-i] . 



So we have the equations 




(22) 



[aki + bki COSTT (x - Xk-i) /Axfc] BkiC^i = 
[ttki-i + bki-i COSTT (x - Xfc_i) /Axfc] Bki-iCl, 
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hence, 

o-kiBkiCli = —aki-iBki-iCli_i, bkiBkiCli = —bki-iBki-iCli_i. 
If we put f l?T]) to these equahties, then we get 

(Zk-ll-l + Zkl-l — Zk-ll — Zkl) f3l = {Zk-ll-l + Zkl-l — Zk-ll-2 — Zkl-2) , 
{Zk-ll-1 + Zkl — Zk-ll — Zkl-l) A = {Zk-ll-1 + Zkl-2 — Zk-ll-2 — Zkl-l) , 

(23) 

where (3i = {Ayi_i/ Ayi) . In analogues the condition of continuity of the 
second derivative F^a on the boundary{xfc_i} x [yi-i, yi] tends to the equations 

[Aki + Bki cosTT (y - yi-i) /Ayi] bkicli = 
= - [Ak-ii + Bk-ii cos 7T{y- yi-i) / Ayi] bk-iicl_^i, 

hence, 

Akibkicli = — Ak-iibk-iicl_ii, Bkibkicli = —Bk-iibk-iicl_ii, 

and 



{Zk-ll-l + Zk-ll — Zkl-l — Zkl) ^k — Zk-ll-1 + Zk-ll — Zk-2l-l — Zk-2h 
{Zk-ll-l + Zkl — Zk-ll — Zkl-l) ^k = Zk-2l + Zk-ll-1 — Zk-2l-l — Zk-ll, 



(24) 



where 6k = (Axk-i/Axk) ■ It is easy to prove that other derivatives of the 
second order are also continuous in this case. 
From dM]) and (jM]) we have 

Zki-2 = l3iZki + Zkl-l (1 - A) , 

Zk-ll-2 = (ilZk-ll + Zk-ll-l (1 - A) , 
Zk-2l = 6kZkl + Zk-ll (1 — 6k) , 
Zk-2l-l = 6kZkl-l + Zk-ll-l (1 — ^fc) • 

Differently these formulas can be written in the form 

{zki-i - Zki-2) IAyf_^ = - {zki - Zkl-l) /Ayf, 
{zk-ii - Zk-2i) = - {zki - Zk-ll) /Axl, 

and this means that 

{zki - Zkl-l) / Ayf = {-l)^ ^.k, l = l,L, k = 0,K 
{zki - Zk-ll) lAxl = i-1)^ ui, k = Tj{, I = Oj: 
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for some sets of values {fik} , {^i}- So, 



I 



Zkl - Zk-ll = Zko - Zk-w + {flk - yUfc-l) ^ (-1)^ A?/| = vi {-if Axl, 



I 



-1) — — 2^(-ir Ay^. = i/,-z/o. 



Axi 

For fulfilment of the last equalities for all fc, / one can take for any value A 



^k 



k 

/ifc = /io + (-ly Ax^l 



In this case 



yi = i/o + A5^(-l)^Ay^ 



2 

i' 



2^00 + ^^0$^ (-1)' Ax2 + fioJ2 (-1)' ^^/|+ 



i=l j=l 

k I 



(25) 



+A5^(-irAa:?^(-iyA2/|. 

i=i j=i 

We see that this family of ^-splines depends only on the next four free pa- 
rameters ZoOjUq, flQ, X. 

Examples of two-dimensional trigonometric splines see in Figs. ITfHl 

4 Superpositions of splines 

In the last section we have established that for two-dimensional sets 971 with 
matrix structure ([2]) solutions of the problems A-C can be found in the form 
of splines (polynomial or trigonometric with different types of smoothness). 
If DJl does not have the matrix form any spline has an additional (false) 
points where its gradient equals to zero; so it can not solve the problem C. 
To overcome this obstacle some constructions may be proposed. 

For example, consider the following procedure consisting of the steps 
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2 4 6 8 

X 

Figure 7: An example of two-dimensional trigonometric spline ( fT9|l . Values 
{zij} , i,j = 1, 10 were taken randomly from a fixed segment. 

Step 1. Complete QJt to a set OJtxy of the matrix form and build any spline 
F (a;, y; SOTxy) that solves the corresponding problem C; 

Step 2. Turn the coordinate system XY for some ("sufficiently small") an- 
gle ipi and build an analogous spline in the coordinates XiYi, 

F {x^,y^;TlxiY^) = F{xcos(fi - y siinpi, x siinpi + y cosLpi;TlxiYi) ] 

Step 3. Repeat the above steps several times for different angles ipk] then 
build the function 

F{x,y) = '^F {xcosLfk - ysmLfk^xsiuLpk + ycosv^fc; OJlx^yJ • (26) 
fc>i 

As gradients of all summands are zero at points 971, the function F {x,y) 
(superposition of splines) solves the problem A. At the same time the addi- 
tional points for matrixes dJlxt.Yk are not the same (they are removed from 
their initial places); so, the probability that gradients of the summary func- 
tion fl26|) at additional points are not zero is high. Test calculations verify 
these assumption. 
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X 



Figure 8: An example of two-dimensional trigonometric C^-spline ( fTOjl . Val- 
ues {zij} , i,j = 1, 10 were calculated by the formula ( [25]) . 

References 

[1] Dolin L.S., Luchinin A.G., Titov V.I., and Turlaev D.G. Correcting im- 
ages of underwater objects distorted by sea surface roughness. - Proc. of 
SPIE Vol.6615. 66150K-1 

[2] Osadchy V.Yu., Savchenko V.V., Levin I.M. , Frantsuzov O.N., and Ry- 
balka N.N. Correction of image distorted by wavy water surface. - Proc. of 
the IV International Conference: Current problems in Optics of Natural 
Waters (ONW'2007), lAP RAS, 2007,pp. 91-93. 

[3] Dolin L.S., Luchinin A.G. Water-scattered signal to compensate for the 
rough sea surface effect on bottom lidar imaging. - Applied Optics, Vol. 
47, No 36, pp. 6871-6878. 



20 



